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ABSTRACT 

We demonstrate that there is an inherent symmetry in the way high-overtone stellar pulsations 
sample the core and the envelope, which can potentially lead to an ambiguity in the astero- 
seismologically derived locations of internal structures. We provide an intuitive example of 
the source of this symmetry by analogy with a vibrating string. For the stellar case, we focus 
on the white dwarf stars, establishing the practical consequences for high-order white dwarf 
pulsations both analytically and numerically. In addition, we verify the effect empirically by 
cross-fitting two different structural models, and we discuss the consequences that this ap- 
proximate symmetry may have for past and present asteroseismological fits of the pulsating 
DBV, GD 358. Finally, we show how the signatures of composition transition zones that are 
brought about by physically distinct processes may be used to help alleviate this potential 
ambiguity in our asteroseismological interpretation of the pulsation frequencies observed in 
white dwarf stars. 
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1 ASTROPHYSICAL CONTEXT 

Asteroseismology is the study of the internal structure of stars 
through their pulsation frequencies. The distribution of the ob- 
served pulsation periods is theoretically determined entirely by the 
run of two fundamental frequencies in stellar models: the buoyancy, 
or Brunt- Vaisala frequency and the acoustic, or Lamb, frequency. 
For nonradial g-mode pulsations in white dwarfs - the only modes 
so far observed in these stars - the buoyancy frequency is the dom- 
inant physical quantity that determines the pulsation periods, re- 
flecting the internal thermal and mechanical structure of the star. A 
considerable body of work by numerous investigators has focussed 
on the signature that core properties, such as the crystallized mass 
fraction and the C/O abundance distributions, and envelope prop- 
erties, such as diffusion profiles and surface layer masses, have on 
the pulsation frequencies. Unambiguous interpretation of the distri- 
bution of periods is hindered by the difficulty of disentangling the 
signatures of structures in the deep core from those in the envelope. 

It is well known both analytically and numerically that the 
separation between the periods of consecutive radial overtones ap- 
proaches a constant value for high-overtone modes, so one might 
think that these periods do not actually contain information about 
the internal structure of the star. The resolution of this apparent 
paradox is that in general there are sharp features in the Brunt- 
Vaisala frequency for which the modes may not be considered to 
be in the asymptotic limit, and the presence of these features per- 
turbs the periods so that the period spacings are no longer uniform, 
an effect which is commonly referred to as 'mode trapping' . Thus, 
these deviations from uniform period spacing contain information 



about sharp features in the Brunt- Vaisala frequency, and can be 
used to discern internal structure in the star, such as the locations 
of composition transition zones. To date, almost all analyses of ob- 
served mode trapping in white dwarf pulsators have focussed on 
determining the thicknesses of the various chemically pure enve- 
lope layers. 

In particular, initial estimat es of the He layer mass of the 
DBV GD 358 yielde d Kr 5 7 M+ feradlev & Wingedll99l . while 
more recent analyses iMetcalfe et all2000ll200lT indicated a glob- 
ally optimal fit with Mh c = 10~ 2 ' 7 M* and a changing C/O profile at 
M, ~ 0.5-0.9M*. However, the Metcalfe et al. analyses still found 
a local m inimum in p ara meter spac e near Mn s ~ 10 _6 M*. Most 
recently, iFontaine & BrassarH J2002) calculated a grid of carbon 
core models that included a double-layered envelope structure, a 
re sult similar to the earlier time-dependent diffusion calculations 
of iDehner & KawaleJ ll995l) . They were able to fit the observed 
periods of GD 358 d own t o a level of precision comparable to the 
fit of IMetcalfe et alj J2001 ). without including any structure in the 
core. Their model had composition gradients at two locations in the 
envelope, near 1CT 3 and 1CT 6 M*, and a chemically uniform core. 

In this paper, we show that there is an inherent symmetry in 
the way in which high-overtone pulsations sample the cores and 
envelopes of the models, with the result that it may be possible to 
reconcile the differences between t he fits mentioned abo ve. Specif- 
ically, by expanding on the work of Montgomery! l2003h . we show 
that a sharp feature in the Brunt- Vaisala frequency in the core of the 
models can produce qualitatively the same mode trapping patterns 
as a corresponding feature in the Brunt- Vaisala frequency in the en- 
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velope, leading to a potential ambiguity in the interpretation of the 
observed periods. In the following sections, we demonstrate how 
this core/envelope mapping works for both DBV and DAV mod- 
els, and we use it to suggest a re-interpretation of the previous and 
current fits for GD 358. As a further example, we show the results 
of cross-fitting between two structurally dissimilar models. Finally, 
we show what considerations allow this 'symmetry' to be broken. 



2 SIMPLE ANALOGY: THE VIBRATING STRING 

In order to better illustrate the problem of g-mode oscillations in 
a white dwarf, we first imagine doing asteroseismology of a far 
simpler object, the vibrating string. As it turns out, all of the major 
results can be carried over to the astrophysical case. 



2.1 The uniform string 

If we take a uniform string with a constant 'sound speed' c, and 
assume a sinusoidal time dependence, then we obtain the following 
eigenvalue equation: 



■ + -^ = 0, <M0) = = </>(L), 



(1) 



dx 2 c 

where xp(x) is the spatial part of the eigenfunction, uo is the eigen- 
frequency, and L is the length of the string. The solution is given 
by 



ip(x) = Asin(k„x), 



HIT 
k„ = T , 



■ k„c, 



1,2,3, 



(2) 



Thus, the spectrum of eigenfrequencies {u} n } is equally spaced. If 
we now make a small position-dependent perturbation to the sound 
speed c, call it Sc(x), then the shift in frequencies may be calculated 
from a variational principle, yielding 



oj„ L 



sin (k n x). 



(3) 



Finally, we may ask what the perturbation to the eigenfrequencies 
is for a new 8c(x) which is the reflection of 8c(x) about the midpoint 
of the string, i.e., Sc(x) = 8c(x), where x = L—x: 



8lu„ 



— I dx 
L 



sin (k„x) 



dx 



{SC(X)\ . 2 

I ——I sin (k„x) 



sin (k„L-k n x) 



(4) 



From comparison of Eqs.|3|and|4| we see that 5c(x) and 8c(x) pro- 
duce identical sets of perturbed frequencies. Thus, from an astero- 
seismological standpoint, the two perturbations are degenerate and 
cannot be distinguished from one another. This is hardly surprising 
since we know that reflection about the midpoint of the string is a 
symmetry of the problem, i.e., it doesn't matter which end of the 
string we call x = 0, and how we choose it had better have no effect 
on the observed frequencies. Thus, the above symmetry must be 
present even if the perturbation 8c(x) is not small. 

2.2 The nonuniform string 

For the case when c = c(x) is a function of x, Eq.Qis still separa- 
ble, but the x-dependence is no longer sinusoidal. However, if we 
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Figure 1. The frequency perturbations due to a bump in the sound speed 
for the nonuniform string: (upper panel) two perturbations which are mirror 
images of one another in the sense given by Eq.|S](solid and dashed curves, 
respectively), both of which produce the same set of frequency perturba- 
tions (lower panel), shown as a function of overtone number n. The vertical 
dotted line in the top panel shows the the geometric location of the 'centre' 
of the string with respect to this reflection. 



consider only modes with relatively high overtone number (n ^> 1), 
then we may use the JWKB approximation to obtain 

%j)„{x) =A—==sin[(/> n (x)], k„(x) = (5) 
Vk„(x) c(x) 

where 

mv J^dx'c' [ (x) 

Wn=-7T , <i>n{x) = n-n— , n= 1,2,3,... (6) 

! a dxc-\x) J Q dx'c-'(x') 

Just as before, a small perturbation to the sound speed, 8c, may be 
related to a change in the eigenfrequencies by 



8lu„ 1 
Un ~ ^dxC-\ X ) Jo 



dxc \x) { ^7^1 sm[(j)„(x)]. (7) 
c(x) ) 



Since the string is not uniform, it no longer possesses reflection 
symmetry about its midpoint. However, in analogy with the uni- 
form case, we consider a reflection coordinate x defined by 



i>„{x) = 0„(L) -</>„(*) 
= im — (f>n(x). 



(8) 



As before, if we consider a 'reflected' perturbation Sc(x)/c(x) = 
8c(x)/c(x), then from Eq.|8|we have dxc~'(x) = -dxc~ l (x), so we 
find that 
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Figure 2. Upper panel: The change in mode periods as a function of over- 
tone number k due to a perturbation to the Brunt-Vaisala frequency at 
M r = 0.5622M* (filled circles) and at a depth of 1(T 5 - 331 M* (open circles). 
Lower panel: residuals of the fit as a function of the position of the envelope 
bump in which the position of the core bump is kept fixed. 



which is identical to Eq.Qif x is replaced by x. Thus, the reflected 
perturbation produces the same set of frequency perturbations as 
did the original perturbation, so there again exists an ambiguity in 
inferring the structure of the string from its eigenfrequencies; this 
is illustrated in Fig.Q 

Finally, we note that this symmetry bears a close resemblance 
to the phenomenon of frequency aliasing in time-series analysis. If 
we consider the frequency perturbation to be a continuous function 
of the overtone number, such as the solid line in the lower panel of 
Fig.0 then the fact that the spectrum of eigenfrequencies is discrete 
means that this curve is sampled only at an evenly spaced set of 
points. Now if the curve in the lower panel of Fig.Qrepresents the 
frequency changes due to a perturbation in 8c /c at x/L = 0.30, then 
the curve corresponding to the symmetric perturbation at x/L « 
0.87 would be a highly oscillatory curve passing through the same 
set of integer values of n. However, since the eigenmodes sample 
this oscillatory curve only at integer values of n, this high frequency 
signal is 'aliased' back to lower frequencies, i.e., the solid curve in 
the lower panel of Fig.Q In this way, a perturbation Sc/c at the 
point x can, through aliasing, appear as if it originates at the point 
X. 



3 WHITE DWARF CORE/ENVELOPE SYMMETRY 
3.1 Analytical approach 

We would now like to apply the ideas of the previous section to 
the pulsating white dwarf stars. We believe this should be pos- 
sible since the adiabatic equations of oscillation of a spherically 
symmetric stellar model (in the Cowling approximation) can be re- 
duced to a second-order differential equation, i.e., the g-mode pul- 
sations can be reduced to a form which mimics that of the vibrat- 
ing string. In particular, the oscillation equation may be written as 
toeubner & Goughll 984 iGoughl 19931) 
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Figure 3. The mapping between points in the core (.r-axis) and those in 
the envelope (v-axis) which produce similar mode trapping for moderate to 
high overtone modes. The solid line is the result of direct numerical fitting 
and the dashed line is the analytical prediction of Eq. 1151 . in which the 
inner and outer turning points, r\ and rj, have been taken to be the center 
and the surface, respectively. 



dr l c 



L - l- N - 

r 2 \ uS 2 



(10) 



where N is the Brunt-Vaisala frequency, L = g(£+l), c is the sound 
speed, and ui c is the acoustic cutoff frequency, which is usually neg- 
ligible except near the stellar surface. For high overtone g-modes, 
we have K ~ LN/ tor, and applying the JWKB approximation we 
find iGouebl 1993b 



and 



kfr) = 



Jfc-i 
2 



f r2 dr'\N\/r''' 



'dr N - 
r 



£=1,2,3,. 



(11) 



(12) 



where r\ and T2 are the inner and outer turning points of the mode, 
respectively, and where we have switched from n to k to denote 
the radial overtone number. The major difference between this case 
and that of the string is that for the string the turning points were 
the same for every mode (i.e., the fixed endpoints), whereas for g- 
modes, the inner and outer turning points of the modes are weak 
functions of the mode frequency, so they are different for every 
mode. Thus, the reflection symmetry may not be as exact as it was 
for the string problem. 

If we consider a perturbation to the Brunt-Vaisala frequency 
8N/N, and we assume that the inner and outer turning points are 
fixed, then the problem is very similar to that of the nonuniform 
string, with the perturbations to the eigenfrequencies given by 

_* r dr m)^\ M r)^]. (B) 
J^drN/rJri \ N J r L 4 J 

Similarly, the approximate reflection mapping between the points r 
and f is given by 



Mr) 



(j>k(ri)-(t>k(r) 



(14) 



which may be explicitly written as 
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dr"- 



\N\ 



r 



(15) 



Essentially, the above equation says that points which are the same 
number of 'wavelengths' (nodes in the radial eigenfunction) from 
the upper and lower turning points are reflections of one another. 

We remark that if we take r\ = and ri = R* in Eq. 1121 then 
<j>k(r) becomes a monotonic increasing function between the centre 
and the surface, and could itself be used as a radial coordinate. In 
particular, if we define 



*(r) = 



JO r 
JO r 



(16) 



then $(0) = and $(_/?*) = 1; since <3> is directly linked to the Brunt- 
Vaisalii frequency, we will call it the 'normalised buoyancy radius' . 
The radial coordinate $ has the advantage that it makes the problem 
look quite similar to that of the uniform string: the 'reflection map- 
ping' between the points r and f is given by <E>(r) = 1 - $(F), with 
the reflection point ('centre of the string') having <£> = |. In addi- 
tion, it can be shown that in the JWKB approximation the kinetic 
energy density per unit $ is constant. Since this approximation is 
valid everywhere except at sharp features such as composition tran- 
sition zones, a change in the kinetic energy density (as well as the 
'weight function') plotted as a function of $ gives us a clear indi- 
cation of the amount of mode trapping which a mode experiences. 
We show examples of this in section|5| 

For completeness, we mention that the analysis of the preced- 
ing paragraphs may also be carried out for high-order p-modes. In 
this case, we have A* ~ uo /c, with $ given by 



fc(r) = 



V dr 

Jr, 



r-dr 

Jr, 



(17) 



Since $ is a function of the sound speed, c, we refer to it as the 
'normalised acoustic radius'. The reflection mapping is unchanged 
from the previous result, namely $(r) = 1 — 5>(r). If a range of I 
values are observed, as is the case for the Sun, then this symmetry 
can be broken due to the fact that the lower turning point, n , is a 
function of £, 



3.2 Full numerical treatment 

In order to investigate this core/envelope symmetry in detail, we use 
a fiducial DB white dwarf model with the parameters M* = 0.6Mq, 
T c ff = 24, 000 K, Mhc = 10~ 6 M*, and we consider pulsation modes 
having I = 1 with periods between 400 and 900 s. We choose the 
diffusion exponents for the C/He transition zone so as to make this 
transition zone as smooth as possible. The reason for this is that we 
wish the background model to be smooth so that the only bumps in 
the Brunt- Vaisala frequency are the ones which we put in by hand. 

In Fig. [5] we show the change in mode periods caused by 
placing a bump in the Brunt- Vaisala frequency at a point in the 
core as well as those produced by a bump in the envelope: the 
filled circles in the upper panel correspond to a core bump at 
M r I Mi, = 0.5622, while the open circles are for a bump in the en- 
velope at log(l -Mr/M^ = -5.331. We see that the shapes of the 
perturbations are qualitatively the same, and that both sets of per- 
turbations are well fit by the same asymptotic formula (solid curve). 
In the lower panel, we show how the residuals change as we vary 
the position of the (induced) envelope bump in order to try to re- 
produce the period changes due to the (intrinsic) core bump, i.e., 
we move the envelope bump through a range of radii in order to 
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Figure 4. Theoretical internal oxygen profiles from Metcalfe et al. 
for a 0.65 Mg white dwarf model produced with standard semiconvec- 
tive mixing (solid) and with complete mixing in the overshooting region 
(dashed) during central helium burning. The two shaded areas show the re- 
gions of the core where perturbations to the Brunt- Vaisala frequency can 
mimic perturbat ions in th e envelope at valu es of log q corresponding to 
those derived bv lFontaine & Brassard l2002t) for GD 358. 



examine whether we have found a global best fit. Clearly, not only 
is there one unambiguous global minimum which is an order of 
magnitude smaller than the other local minima, the residuals are a 
smooth function of the bump position, so we are justified in using 
a nonlinear fitting algorithm in order to find the local (and global) 
minimum at log(l -M r /M+) = -5.331. 

Using the above procedure, we can map out the corresponding 
pairs of core/envelope points in our equilibrium model, obtaining 
the solid curve in Fig. [3] (the dashed line is the result of using the 
analytical relation given by Eq. d 1 5i with r\ = and ri = /?*). Thus, 
we see that a feature in the core at M, -/M* ~ 0.5 can mimic a feature 
in the envelope at log(l —M r /M t ) ~ -5.5. This is a very significant 
result because we a priori expect there to be both envelope bumps 
due to chemical diffusion as well as core bumps due to the prior 
nuclear burning history of the white dwarf progenitor. 

For instance, in Fig. [4] we show evolutionary C/O profiles 
ISalaris et alJll997l) which result from assuming either standard 
semiconvective mixing (solid line) or complete mixing in the over- 
shoot region (dashed line). Since the Brunt- Vaisala frequency de- 
pends on the radial derivative of these profiles, the regions of large 
slope at M r /M„ ~ 0.5 and M r /M„ ~ 0.98 will produce bumps in 
the Brunt- Vaisala frequency at these points. For the shaded regions 
in this figure, we have taken the quoted ran ges of the He t ransi- 
tion z ones in the envelope of the model of iFontaine & Brassard 
(2002) and used the reflection mapping of Fig.|3|to indicate which 
regions in the core correspond to these envelope ranges. Almost 
eerily, these ranges correspond quite closely to those in which we 
would e xpect to see str ucture in the core. Thus, the possibility ex- 
ists that lFontaine & Brassard are fitting 'real' structure in the core 
with assumed structure in the envelope. The reverse is also possi- 
ble, of course, and at the very least the potential exists for the two 
signatures to become entangled with one another. 
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Figure 5. The same as Fig.|2| but for a DAV model. From the bottom panel, 
we see that the local minimum near ~ —8.2 (arrow 2) has residuals nearly 
as small as those at the minimum near ~ —7.2 (arrow 1), and that the resid- 
uals near ~ -6.4 are also fairly small. Thus, for all practical purposes, the 
core/envelope mapping is no longer single valued. 

3.3 DA models 

The above analysis can be directly applied to DAV models. For 
our fiducial model, we take M* = 0.60M Q , % s = 12, 000 K, M He = 
10~ 3 M*, Mh = 10"* M*. In order to make our background model as 
smooth as possible, we ignore the Ledoux term in the computation 
of the Brunt- Vaisala frequency; this reduces but does not eliminate 
the bumps due to the chemical transition zones. Finally, since these 
stars are cooler and therefore more degenerate, we consider longer 
period 1=1 modes, of between 800 and 1300 s. 

For the DAVs, we find a couple of surprises. First, as shown 
in Fig. [5] the mapping is no longer single-valued: there is more 
than one point in the envelope which is approximately mapped to a 
point in the core, and vice versa. In fact, there appear to be at least 
two and possibly more families of mappings which have approx- 
imately the same residuals. This is due to the fact that our back- 
ground model already has two 'bumps' in it corresponding to the 
C/He and He/H transition zones, which our test bump is interacting 
with. 

Second, as shown in Fig. |6| the core/envelope mapping is 
weighted much more toward the envelope than it is for the DBV 
models, i.e., a given point in the core is mapped to a point farther 
out in the envelope than it is in the DBV models. Physically, this is 
because the DAV model is cooler and has a more degenerate core, 
and therefore N is smaller in its core. Thus, one has to move far- 
ther out from the centre in order to accumulate a given amount of 
phase in the sense of Eq. I15i or <16> , The implications of this could 
be important, since a changing C/O profile in the range 0.5-0.9 M t 
maps to an envelope point of 10~ 10 -10~ 8 . Thus, such a C/O profile 
cannot mimic a hydrogen layer mass greater than ~ 10~ 6 M*, so the 
ambiguity in interpreting the origin of perturbations to the periods 
may be alleviated in the DAVs. 

Finally, we wish to mention that at least one of the DAVs, 
BPM 37093, is theoretically predicted to have a partially crystal- 
lized core. In terms of seismology, the main effect of such a core 
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Figure 6. The same as Fig. [3] but for a DAV model. The dashed curve is the 
analytical result for the fiducial DAV model, and the two solid lines labelled 
1 and 2 are the core/envelope mappings corresponding to minima 1 and 2 in 
Fie.pl For reference, the dotted curve is the analytical result for our fiducial 
DBV model. 

would be to exclude the pu lsations from the crystallized region 
iMontgomery & WingeJ [l999albl) . effectively making the lower 
turning point, n, the upper boundary of the crystallized core. Thus, 
instead of running from to 1, the horizontal axis in Fig. |5| would 
run from Afe/M*, to 1 , where M&/M+ is the crystallized mass frac- 
tion. This would tend to push the symmetry point for a given loca- 
tion in the core farther out into the envelope. Also, the higher mass 
of BPM means that its core will be more degenerate than that of the 
other DAVs, which will also push its symmetry mapping yet farther 
out into the envelope. 



4 CROSS-FITTING TWO STRUCTURAL MODELS 

Empirical evidence that this symmetry in the white dwarf mod- 
els may lead to some ambiguity in the interpretation of model-fits 
from var ious descriptions of the stel lar in terior was recen tly pub- 
lished in lFontaine & Brassard 120021) and lMetcalfel <2003l) . To in- 
vestigate the effects of this symmetry more directly, we performed 
several cross-fitting experiments with two different structural mod- 
els - attempting to match the calculated pulsation periods from one 
model by using a structurally distinct model to do the fitting. In 
particular, we attempted to fit the carbon core double -layered en- 
velope model periods from Table 1 of Fontaine & Brassard 1 2002) 
using the 5-parameter model of Metcalfe et al. 1 2001 ), which in- 
cludes a single-layered envelope and an adjustable C/O core; we 
were able to find a match with root-mean-square period residuals of 
only <tp = 1.26 second s (the opt i mal m odel parameters were iden- 
tical to those found bv lMetcalfd <2003l) for GD 358). 

We may also use Fontaine & Brassard's published results 
to turn the question around and see what goodness of fit they 
would obtain by attempting to fit the periods of our model. An 
upper limit on the root-mean-square residuals that may arise from 
such a comparison can be obtained by comparing the perio ds of 
the two optimal models fo r GD 358 from Metcalfe I2003J) and 
iFontaine & Brassarcl 120021) . leading to ap ^ 1.26 seconds. It is 
important to note that the double-layered envelope models may be 
able to match the periods from the adjustable C/O models even bet- 
ter than this with slightly different model parameters - a possibility 
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that would even more clearly demonstrate a core/envelope symme- 
try. 

Given the discussion of the previous sections, the results from 
the cross-fitting are not too surprising , since we are in e ssence fit- 
ting structure in the envelope of the Fontaine & Brassard model us- 
ing structure in the core of our model located at the appropriate 
reflection point, and vice versa. The level of the residuals from 
our model-to-model comparison tells us directly what we have 
suspected for some time: the model-to-observation residuals for 
GD 358 (up ~ 1 second) are dominated by structural uncertain- 
ties in the current generation of models, regardless of which type 
of model is used to do the fitting. This does not necessarily mean 
that the conclusions based on fitting from either of these models 
should be thrown out; it simply means that neither model is a com- 
plete description of the actual white dwarf stars, a statement that 
can hardly be considered controversial. 

We can attempt to determine which of the two structural 
descriptions is closer to reality (or better describes the interior 
structure as sampled by the pulsations) by comparing the abso- 
lute level of the residuals for GD 358 from each model, cor- 
rected for the nu mber of free parameters. The 4-parameter model of 
iFontaine & Brassarcl 12002) leads to op = 1.30 seconds when com- 
pared to the periods observed in GD 358. The Bayes Information 
Criterion jKoen&Lanevfe oOO) would lead us to expect the residu- 
als of a 5-parameter fit to be reduced to 1.17 seconds just from the 
additio n of an extra parameter. The 5-parameter model of Metcalfe] 
l2003h does much better than this, with residuals of op = 1.05 sec- 
onds (an improvement equivalent to 4a b s ) - suggesting that the 
internal C/O profile may be the more important of the two possi- 
ble structures. Based on the results of our own experiments with 
double-layered envelope models (Metcalfe et al., in preparation) 
this possibility may be even more likely. 



5 'BREAKING' THE SYMMETRY 

In the previous sections, we have demonstrated that a core/envelope 
symmetry exists for high-overtone modes. Fortunately, this sym- 
metry is approximate, and there are many ways in which it may 
be lifted or broken. For example, many of the DBVs appear to be 
higher overtone pulsators (k ~ 10), although this is not necessar- 
ily true for all of the observed modes in a given star, nor is it true 
of every member of the class. The DAVs as a class are not high- 
overtone pulsators, so this symmetry will be less of an issue for 
them. In addition, as shown in Fig. |6| the core/envelope mapping 
for DAV models is weighted more towards the envelope, so that the 
reflection point for many points in the core is in a region of the en- 
velope where it may not overlap with the expected position of the 
chemical transition zones. 

In addition, the core/envelope symmetry can be broken if we 
make additional assumptions. For instance, due to the different 
physical processes which produce them, the generic shape expected 
for the C/O profile in the core will be different from the expected 
shape of the C/He profile in the envelope. Using this information 
(parametrized in some form, for instance), we should be able to dis- 
cern a core feature from an envelope feature. We show an example 
of this in Fig. Q in which we see that the reflected bumps (lower 
panel) do indeed have different amplitudes and shapes from the 
actual bumps (middle panel) in the Brunt-Vaisala frequency [note 
that we have used the buoyancy radius, $, as defined in Eq. 1161 as 
our radial coordinate; along the top axis we use the more familiar 
log(l-M r /M*)]. 
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Figure 7. The bumps in the Brunt- Vaisala frequency (middle panel) which 
are produced by given chemical transition zones (upper panel), and the mir- 
ror image of these bumps under the reflection mapping (lower panel). We 
have taken the buoyancy radius <& to be the radial variable; along the top axis 
of the upper panel we indicate the corresponding values of log(l — M r /M*). 



The physical inp uts used to generate Fig.0include a Salaris- 
like C/O profile (e.g., Salaris et al. 1997) and a two-tie red C/He 
Dehner-like envelope profile l Dehner & Kawalen ll995l) . Further- 
more, we have defined a 'bump', SN/N, as the fractional differ- 
ence between the Brunt-Vaisala frequency calculated both with and 
without the Ledoux term (this term explicitly takes account of the 
effect which composition changes have on the value of N ). While 
this yields reasonable results for the inner transition zones, for the 
outer C/He transition zone near log(l -Mr/M*) ~ -6 this prescrip- 
tion would greatly overestimate the importance of this bump. In- 
stead, for this case we have taken SN/N to be the fractional dif- 
ference of the actual Brunt- Vaisala frequency and a third-order 
polynomial which smoothly joins onto the Brunt-Vaisala frequency 
on each side of the transition zone. With these definitions, we see 
that the inner C/O transition zone should produce the most domi- 
nant mode trapping feature, and that the inner C/He transition zone 
should produce the next most important feature, since these bumps 
are both the highest and the narrowest features present. In contrast, 
the outer C/He transition zone should have a relatively minor effect 
on the mode trapping. 

We illustrate the mode trapping ability of this model in Fig.[S] 
in which we plo t, again as a fu nction of $, the weight functions 
for/V 2 (see lKawaler etalJll985l their Eq. 8c) of the full numeri- 
cal problem: the middle panel is the weight function for an I = 1, 
k = 1 1 mode, and the lower panel is that of an ^=1, £ = 21 mode. In 
the upper panel we show the bumps corresponding to the composi- 
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Figure 8. Examples of mode trapping of the eigenfunctions due to the com- 
position transition zones. The upper panel shows the bumps in the Brunt- 
Vaisala frequency from Fig. Q as a function of the buoyancy radius, 
and the lower two panels show the weight function WF§ for two different 
modes, k = 1 1 and k = 2i (both with t = 1). The shaded regions indicate the 
extent of the two most important transition zones; the differing amplitudes 
on each side of the transition zones clearly illustrates the mode trapping 
which these zones produce. 

tion transition zones, and in all panels we have shaded the regions 
containing the two largest bumps, since these should produce the 
strongest mode trapping. As we can see, this is certainly the case 
for the two modes we have selected: the modes show a marked 
change in amplitude in these shaded transition regions, but other- 
wise they propagate with essentially 'constant amplitude', i.e., their 
amplitudes evolve according to the JWKB formula. 

So how does mode trapping help break the core/envelope sym- 
metry? First, imagine that we wish to calculate the shift in peri- 
ods due to all the bumps, not just to the two large ones we have 
highlighted. The k = 21 mode is trapped between the two transi- 
tion zones, so it will sample the bump at $ ~ 0.37 more strongly 
than modes such as the k = 11 mode, which is not trapped in this 
region. Conversely, the k = 1 1 mode has an enhanced amplitude in 
the outer layers ($ > 0.6), so it will sample the bump at $ ~ 0.77 
more strongly than modes such as the k = 21 mode, which is not 
trapped in this outer region. Thus, the effect which the smaller 
bumps have on the frequencies is strongly influenced by the mode 
trapping which the larger bumps produce. As a result, it should in 
principle be possible to discern whether there are in fact two bumps 
rather than one large one, what their relative locations are (e.g., both 
in the core, both in the envelope, or one each in the core and enve- 
lope), and what their relative strengths are. 



Additionally, the presence of modes of different I may help 
in resolving this core/envelope symmetry, since the outer turning 
point is a function of I for many of the mod es. Also, const raints on 
the rotational splitting kernels (e.g.. iKawaler et al]ll999h derived 
from rotationally split multiplets may be of help. In future calcula- 
tions, we will attempt to address quantitatively the conditions un- 
der which it is possible to resolve the core/envelope symmetry, both 
in terms of the number of modes required as well as the physical 
assumptions needed regarding the shapes of the features to be re- 
solved in the Brunt-Vaisala frequency. 



6 DISCUSSION & CONCLUSIONS 

We have demonstrated, both analytically and numerically, that a 
symmetry exists which connects points in the envelope of our mod- 
els with points in the core of our models. Specifically, we find that 
a sharp feature ('bump') in the Brunt-Vaisala frequency in the en- 
velope of our models can produce the same period changes as a 
bump placed in the core, and we have numerically calculated this 
core/envelope mapping. While we have restricted ourself to the 
case of high-overtone white dwarfs which pulsate in g-modes, such 
a symmetry should be a generic feature of all high-overtone pul- 
sators, whether they pulsate with g- or p-modes. 

The specific motivation for much of our analysis has been the 
well-studied DBV, GD 358. Given its mass, stellar evolution theory 
leads us to expect that it has a C/O core, and that in some region of 
the core there should be a transition from a C/O mixture to a nearly 
pure C composition. If this transition begins at ~ 0.5M*, then it 
will produce a bump in the Brunt-Vaisala frequency at this point. 
As discussed in the previous sections, such a bump could be mim- 
icked by an envelope transition zone with a depth of ~ 10~ S M*. 
Therefore, any incompleteness in the modelling of the core bump 
could be 'corrected' by a transition zone in the envelope placed at 
a depth of 10~*'°M*. This could explain the initial He layer deter- 
mination of lff^M* by Bradley & Winget (1994) as well as the 
continuing presence of a loc al minimum near ~ MT 6 ' M* in the 
current generation of models I Metcal fe" et al 1200(1.120011) . 

In addition, there is reason to believe that the actual C/He 
profile may be two-tiered, with a transition from pure C to 
a C/He mixture at ~ 10~ 2 M* and another transition to pure 
He at ~ 10~ 6 M+ iPehner & KawaleJll995t ICorsico et all l2002t 
iFontaine & Brassard^OO^ : if this is the case, it is quite possible 
that the mode trapping effects of a C/O transition zone in the core 
could become entangled with those of the outer He transition zone. 
While additional considerations may break this core/envelope sym- 
metry, we need to be aware of this aspect of the problem in order 
to make progress in modelling these stars. 

In summary, we have shown that, for moderate to high over- 
tone pulsators, there exists a symmetry in the mode trapping pro- 
duced by features in the core and those in the envelope which can 
lead to ambiguity in determining the location of features such as 
composition transition zones. This may explain the present and pre- 
vious fits for GD 358, and at the very least is something which must 
be taken into account in future asteroseismological fits. 
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